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Abstract 

In the context of massive gravity, bi-gravity and multi-gravity non-minimal 
matter couplings via a specific composite effective metric were investigated 
recently. Even if these couplings generically reintroduce the Boulware-Deser 
ghost, this composite metric is unique in the sense that the ghost reemerges 
only beyond the decoupling limit and the matter quantum loop corrections 
do not detune the potential interactions. We consider non-minimal derivative 
couplings of the composite metric to matter fields for a specific subclass of 
Horndeski scalar-tensor interactions. We first explore these couplings in the 
mini-superspace and investigate in which scenario the ghost remains absent. 
We further study these non-minimal derivative couplings in the decoupling- 
limit of the theory and show that the equation of motion for the helicity-0 
mode remains second order in derivatives. Finally, we discuss preliminary 
implications for cosmology. 


1 Introduction 

The accelerated expansion of the late Universe is confirmed by many cosmological 
observations like supernovae, CMB, baryon acoustic oscillations and lensing. It is 
believed that this expansion is driven by dark energy, whose origin has not yet been 
identified. There are many attempts to explain the nature of dark energy. Assuming 
that General Relativity is the ultimate theory for gravitational interactions, one can 
relate the effective repulsive force between cosmological objects with the presence of 
a cosmological constant A with a constant energy density. In fact, from an effective 
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field theory point of view it must be included. From a particle physics perspective 
the origin of the cosmological constant could be related to the vacuum energy den¬ 
sity. Nevertheless, the theoretical expectations for it using standard techniques of 
quantum held theory exceeds the observational bounds by a tremendous amount. 
This constitutes the century problem of modern cosmology [1, 2]. 

As an alternative explanation for the recent acceleration of the Universe one 
can invoke new dynamical degrees of freedom, either in form of new huids with 
negative pressure in the stress energy tensor, or in form of modihcations of the 
geometrical part of Einstein’s equations. The main goal of these approaches is to 
weaken gravity on cosmological scales. These new degrees of freedom could play 
the role of a condensate, whose energy density either sources self-acceleration or 
compensates the cosmological constant. Such scenarios could for instance arise in 
massive gravity or in higher-dimensional frameworks. 

Concerning the higher dimensional framework, the Dvali-Gabadadze-Porrati (DGP) 
model was a promising candidate among large scale modihed theories of gravity 
[3]. In this scenario our Universe is conhned to a three-brane embedded in a hve- 
dimensional bulk. The brane curvature in form of an intrinsic Einstein Hilbert term 
guarantees the recovery of General Relativity on small scales. On the other hand, 
on cosmological scales the modihcations can be interpreted as a soft mass of the 
graviton, that systematically limits the effective range of the graviton. A mass¬ 
less higher dimensional graviton carries the same number of degrees of freedom as 
a massive graviton on four dimensions, hence the effective graviton on the brane 
carries the usual helicity-2 modes, two helicity-1 modes and one helicity-O mode as 
in massive gravity. The phenomenological consequences of these theories is tightly 
related to the helicity-O mode, which can mediate an extra hfth force. The helicity-1 
modes usually decouple. Even in the presence of this helicity-O mode one can re¬ 
cover General Relativity via the Vainshtein mechanism [4]. It relies on the successful 
implementation of screening the unwanted modes from the gravitational dynamics 
on small scales via nonlinear interactions. The coupling of the helicity-O mode to 
matter is suppressed once one canonically normalises the helds. This decoupling of 
the nonlinear helicity-O mode is manifest in the so-called decoupling limit. Thanks 
to this helicity-O mode, the DGP model admits self-accelerating solutions, which 
however are plagued by ghost-like instabilities [5, 6]. This gave rise to attempts for 
considering more exhaustive setups with Gauss-Bonnet terms in the bulk [7]. 

The decoupling limit of DGP contains a very specihc derivative interactions for 
the helicity-O mode. Their form is dictated by the residuals of the 5-dimensional 
Poincare invariance. Soon it was realised that the decoupling limit of DGP can be 
extended to more general interactions with invariance under internal Galilean and 
shift transformations [8]. Together with the postulate of second order equations of 
motion these symmetries severely restrict the allowed derivative interactions. These 
Galilean invariant interactions have been treated in their own right as scalar held 
theories and received much attention for their rich phenomenology [9, 10, 11, 12, 
13]. In the presence of a compact source these Galileon theories admit spherically 
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symmetric solutions that propagate superluminal [14, 8, 15], even though a direct 
analogue of Hawking’s chronology protection conjecture prevents to construct closed 
timelike curves [16]. The generalisation of the Galileon to a non-flat background 
breaks the Galilean symmetry. It was also realised that one has to invoke non- 
minimal coupling between the scalar held and the curvature in order to maintain the 
property of second order equations of motion [17]. The same covariant interactions 
can be also obtained from a higher-dimensional construction [18]. Furthermore, 
the Galilean symmetry can be promoted to a generalized Galileon symmetry on 
maximally symmetric backgrounds. For instance, the de Sitter Galileon interactions 
acquire additional, specihc potential-like terms [19, 20]. 

Interestingly, these Galilean invariant interactions naturally arise in gravitational 
theories with a massive graviton as an exchange particle [21, 22]. The existence of 
a graviton mass is a fundamental question from a theoretical perspective. It was 
a challenge to construct a consistent non linear theory for massive gravity without 
the Boulware-Deser ghost, which has witnessed a tremendous amount of progress 
in the past decade [21, 23, 24, 25, 26, 27], and shown to be technically natural 
[28, 29]. Besides the interesting analogue to Galileon interactions, one can also 
connect massive gravity to Horndeski type of theories. Horndeski interactions are 
the most general scalar-tensor interactions which give rise to second order equations 
of motion [30, 31]. The covariantization of the decoupling limit of massive gravity 
yields a proxy theory to massive gravity [32], which is a non-minimally coupled 
subclass of Horndeski scalar-tensor theory, which contains four arbitrary functions 
in the Lagrangian. In the proxy theory, these four arbitrary functions are dictated by 
the covariantization procedure, and it shares the same decoupling limit with massive 
gravity. The phenomenological implications of this proxy theory were studied in 
[32, 33]. 

The Galileon interactions are non-renormalizable. Exactly the same property 
is also fulhlled by the decoupling limit interactions of massive gravity [28, 34]. Of 
course this is not the case anymore when one considers also matter loop quan¬ 
tum corrections. In this context, it is crucial to investigate the possible consis¬ 
tent couplings to matter helds in massive gravity. One important requirement is 
to couple the matter fields in a way that does not detune the potential interac¬ 
tions at the quantum level. There is an unique composite effective metric built 
out of the two metrics that fulhls this requirement and guarantees the absence of 
the Boulware-Deser ghost in the decoupling limit [35, 36, 37]. The presence of this 
effective composite metric has important theoretical and phenomenological conse¬ 
quences [35, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55]. In 
this work we would like to further study the theoretical implications of this effec¬ 
tive composite metric in the context of non-minimally coupled scalar-tensor theories 
motivated by the proxy theory. For this purpose, we consider the specihc scalar- 
tensor interactions of the covariantized decoupling limit, that live on the effective 
composite metric. It is important to investigate whether non-minimal couplings to 
this effective metric maintain the ghost absent in the decoupling limit, and whether 
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this might yield richer cosmological phenomenology than the standard Horndeski 
type of theories living in the standard metric. 

2 Non-minimal derivative couplings 

An interesting question within massive (hi-) gravity that has received much attention 
lately is the consistent coupling to matter helds. In this context a new composite 
effective metric was proposed in [35] and it was shown, that it is the unique metric 
that removes the Boulware-Deser ghost in the decoupling limit and through which 
the matter loop quantum corrections do not detune the potential interactions. For 
an application to dark matter, in [54] it was shown that the kinetic Lagrangian 
needs to be restricted to either Mg^—gRg + as in the standard ghost- 

free kinetic terms in bigravity with dynamical g and / metrics [56], or alternatively 
one can replace one of these kinetic terms by the corresponding kinetic term of 
this new composite effective metric M‘^y/—gRg -|- —ges_Rea- It was proven 

that one has to restrict the theory to have not more than two kinetic terms in 
order not to reintroduce the BD ghost in the decoupling limit. In this work, we 
go along the line of investigations as in [54], but aim to push forward this analysis 
to non-minimal derivative couplings to matter helds and investigate whether these 
derivative couplings reintroduce the ghost in the decoupling limit and weather the 
ghost propagates on top of maximally symmetric backgrounds. For this purpose we 
choose to work with the following action 

“S' = y y/^Mp^R + -h y/-gesC^iges, tt) . (1) 

with non-minimal derivative couplings in the matter Lagrangian C.^ and the poten¬ 
tial interactions in form of y/ges- The composite effective metric has a very specihc 
structure in terms of the g and / metrics [35] 

9‘iH = + 2a/3 , (2) 

with a and /3 being two arbitrary constants, and stands for 

Of. = 9 ,pK = f,pYP , (3) 

where we denoted the square root matrix by X = \/9~^f, and its inverse by V = 
\/ f~^ 9 . The special property of this effective metric is the fact that its square root 
of the determinant corresponds to the allowed potentials 

\/-9eS = det(al -f- (3X) = \fOYf det(/91 aV) . (4) 

We can express the square root of the determinant also in terms of the elementary 
symmetric polynomials en(X) and e„(X) as 

4 4 

n=0 n=0 
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where the elementary symmetric polynomials have the following specihc structure 


eo(A') = 

1. 

(6a) 

ei(A) = 

W. 

(6b) 

e,(X) = 

i(W^-[Ay). 

(6c) 

e,(X) = 

i([A]^-3[A][AT+2[A’]). 

(6d) 

e,(X) = 

i ([A] ‘ - 6 [A] V 3 [A=] ^ + 8 [A] [A V 6 [A^) . 

(6e) 


The Lagrangian for the non-minimal derivative couplings between the composite 
effective metric and the scalar held is given by 

C, = A4, . (7) 

In Lagrangian (7), is the Einstein tensor of the effective metric and the tensor 
corresponds to the dual Riemann tensor of the effective metric 



Coupling the matter held to this ehective metric (instead of coupling it only to the 
space-time metric g) has important theoretical and phenomenological consequences. 
Particularly, it is a very interesting question to study the consequences of non- 
minimal derivative couplings to this ehective metric. We could in principal consider 
all the Horndeski scalar-tensor interactions. However, in this work we would like 
to focus on a subclass of the Horndeski interactions given by (7). This subclass is 
special in the sense that it is naturally generated by the covariantisation procedure 
of the decoupling limit of massive gravity discussed in [32, 33]. Note, that we did not 
include an additional kinetic term for the / metric in our Lagrangian (1) in contrast 
to the standard formulation of bigravity [56]. The reasons for that will become 
clear in the next section, when we perform the Hamiltonian analysis of maximally 
symmetric space-times in the context of allowed kinetic terms in the presence of 
derivative non-minimal couplings to the matter helds. 

As next, we shall compute the covariant equations of motion. First of all the 
Einstein equations are given by [25, 43]^ 


2Af|, rye,"'" - ^ 

n=0 


n=0 


^We will follow the same notation as in [54]. 
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where the tensors and stand for [25] 


U^n) = Y.^-r^n-r>{X)Xr 

p=0 
n 

v,n) = ^(-)’’e„_,(y)y^ 

p=0 

and the stress energy tensor is given by 

T^Seff _ ^ J^SeB(2) _ _^y^3eff(3) 

pu pu ^3 fJ-i' ^6 

Note that the properties of the Einstein and Riemann dnal tensor ensnre that the 
scalar held vr enters at most with two derivatives in the stress-energy tensor, 

7 -' 3 efr(l) _ ^ 9 eff(l) TT/gfSeff 

flU ' flU ' 

Q 

7^9eff(3) _ 9efr(3) j T9eS !l9i 

~ + 2 ^ ^ 

with X = (c?7r)^ and the tensors encode the interactions of the scalar held 


^9esW 

= |n«]9;:-n“ 

(13) 

.7-9eff(2) 

= (n“)" - |n,.,]n« - i(|nj,j - [Yi,„Y)g^ 

(14) 

.7-9eff(3) 

= 6(n«)“ - 6in..,|(n«)= + adn,.,]^ - |nj_j)n« 



- 3|njj|n„] + 2iny), 

(16) 


with the abbreviation [ff^^jj] = for convenience. The eqnation of motion 

for the scalar held vr is given by 

£. = (2n«n« + = o. (i6) 

where we made nse of 

vrvrvr’^ir.'"’ = (v) 

with the Ganss-Bonnet term of the ehective metric £q§ = -|- {R^^au/sY ~ 

In the following we wonld like to stndy the conseqnences of the pres¬ 
ence of the derivative non-minimal conplings in the action (1). 


(10a) 

(10b) 

( 11 ) 
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3 Hamiltonian analysis 


In this section we would like to perform the Hamiltonian analysis in the mini¬ 
superspace and discuss different scenarios. For this purpose we shall concentrate on 
the first non-minimal coupling term 'n'Rg^g in (7). The conclusions we make here 
for this first interaction also apply to the remaining interactions. We will consider 
three different scenarios concerning the dynamics of the / metric appearing in g^fi. 

1 . Non-dynamical / metric: The hrst scenario corresponds to the case where 
the / metric is not dynamical and we choose it to be Minkowski. Thus, we have 

ds^ = —N{tYdt‘^ -I- a{tYdx‘^, and dsj = —dt^ + djp . (18) 

As we mentioned above, we will focus on the sub Lagrangian 

£(1) = Mp^y/^Rg - MeffOi . (19) 


but all the statements that we make here apply also to the remaining non-minimal 
couplings of our Lagrangian (7). The conjugate momenta associated with the tt field 
and the scale factor are 


Ptt = 


GMesaiaalffO, 


N, 


eff 


Pa 


12Mp^ad 

N 


+ 


QaiMesaaes{2a7rd -|- OeffTr) 


N, 


eff 


( 20 ) 


with the effective lapse function N^sit) = aN(t) -|- {3 and the scale factor aefr(t) = 
aa{t) + (3- After performing the Legendre transformation we obtain that the Hamil¬ 
tonian has a piece that is non-linear in the lapse 

This is a clear indication that we have lost the constraint equation that removes the 
Boulware-Deser ghost. The presence of the derivative non-minimal coupling with 
non-dynamical / metric yields an unavoidable ghost degree of freedom. This is in 
complete agreement with the findings in [46], where they considered couplings in 
form of Horndeski scalar-tensor interactions living in the composite effective metric. 
In the case they studied, the / metric was not dynamical and their perturbation 
analysis about a FLRW background indicated the presence of the ghost degree of 
freedom. Here we have seen that the Hamiltonian is always non-linear in the lapse 
in the presence of this type of derivative non-minimal couplings. However, as we 
shall see as next, we gain more freedom in making the / metric dynamical as well 
and the Hamiltonian becomes automatically linear in the lapses if we do not add an 
additional kinetic term for the / metric. Actually, this same conclusion applies to 
the pure kinetic Lagrangian without the non-minimal derivative coupling. If one for 
instance takes the kinetic terms M^^—gRg -|- ^fegReff with non-dynamical / 

metric, then this restricts the phase space so much that the Hamiltonian becomes 
non-linear in the lapse. If one allows the / metric to be dynamical but without 


( 21 ) 
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including its kinetic term, this enriches the phase space and the Hamiltonian has 
only then the chance to be linear in the lapses [54]. 

2. Dynamical / metric without additional kinetic term: In this scenario 
we assume that the / metric also carries dynamics in but without including an 
additional kinetic term for the / metric. This is a crucial distinction between this 
scenario and the next one. So our starting point is again the sub Lagrangian (19) 
but this time with 

dSg = —N{tYdt^ + a{tYdx^ , and ds'j = —+ a/(t)^dx^ . (22) 

In this case the conjugate momenta associated with the tt held and the two scale 
factors become 


eaiMeffa^jjaeff UM^^aa 

Pn = - p^^ - 


iVeff 

where we have dehned 


N 


Q = 


6ci]^TfgffCigff (27r(X0ff ~\~ n0ff7r) 


N, 


eff 


This time the Hamiltonian 

= dpg + dfTif + 7rp^ — , 


(23) 


(24) 


(25) 


is automatically linear in the lapses. Thus, giving dynamics to the / metric res¬ 
cues the derivative non-minimal couplings. The linearity of the Hamiltonian in the 
lapses ensures the propagation of the constraint that on the other hand removes the 
Boulware-Deser ghost. 

3. Dynamical / metric with an additional kinetic term: In this third 
scenario we shall consider the dynamical case as in the second scenario but with an 
additional kinetic term for the / metric. This has severe consequences. Our sub 
Lagrangian this time consists of 

£(3) ^ . (26) 

In this scenario the conjugate momenta for the scale factor of the / metric obtains 
an additional term while the others remain the same 


Ptt = 


6aiM0fjagfja0ff 




eflf 


Pf 


12Mjafdf 


+ /5Q, Pg = 


12M, 


PI 


aa 


N 


+ aQ (27) 


After expressing the time derivatives of the helds in terms of their conjugate mo¬ 
menta, one obtains that the Hamiltonian is always non-linear in the lapses inde¬ 
pendently of what one choses for the coefficient in front of the non-minimal cou¬ 
pling. This non-linear contributions are proportional to In order for 








the Hamiltonian to remain linear in the lapses in the presence of the derivative 
non-minimal couplings, one has to restrict himself to the second scenario. It is 
mandatory to have both metrics dynamical but without the inclusion of the second 
kinetic term. This coincides with the hndings of [54] where it was shown that one 
can not have three kinetic terms at the same time. Here we saw also that exactly the 
same happens in the case of derivative non-minimal couplings to matter helds. We 
have checked explicitly that these conclusions also apply to the remaining derivative 
couplings y/-gesdf,7rd^7TGi^^^ and . 

Summarising, among the three different scenarios that we have considered here, 
only the second option yields a consistent way of coupling. If one insists on non- 
minimal derivative couplings of the composite effective metric to the matter helds, 
then this excludes the more restrictive case of hxed / metric. It would be very 
interesting to perform the same cosmological perturbation analysis as in [46] but 
keeping the second metric dynamical as well without the inclusion of its kinetic 
term. 


4 Decoupling Limit 

In this section we will pay special attention to the non-minimal couplings between 
the scalar held vr and the ehective metric in the decoupling limit. First we restore 
the broken diheomorphism by introducing the Stueckelberg helds in the / metric 

f/ii' f/iu = fabdiiS°‘dyS^. (28) 

These Stueckelberg helds 8°" can be further decomposed into the helicity-0 mode 0 
and helicity-l mode 


mA°‘ 

aT 


(29) 


Now we can take the decoupling limit consisting of sending Mpi —)■ oo, while keeping 
the scale Ag = Mpim^ hxed. Let us mention that we will neglect the contribution 
of the helicity-l held and keep track of the contributions of the helicity-0 mode 0 
only. We would like to see whether or not the equation of motion for the helicity-O 
mode contains higher derivative terms due to the derivative non-minimal couplings. 
For this, we will concentrate on the part of the action (1) and we will further 
assume in this section and 

^ fIJ.V i.VfJ.l' y (30) 

where stands for In this case the ehective metric in the 

decoupling limit becomes 


afu ~9fv = [(« + ^ • 


(31) 
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The contribution of the derivative non-minimal couplings C.,^ to the equation of 
motion for the helicity-0 mode is given by 


5C^ _ f 6C^ 5~gf^ \ 

5(t> A3 

= -% {m\nuTX) [(a+ m - f^^'A - (32) 

The stress energy tensor is tranverse on-shell 

(\/^r£) + rj|.« rj;; = . ( 33 ) 

where 8-,^ denotes the equation of motion of the matter held tt given by (16). Using 
this we can write the equation of motion for the helicity-0 mode as 

^ {/3 (\/=l^y>r£. - \/=i;^r;rTr) [(a + P)Sl - /St-] 

This can be similarly written as 

where 

T^‘'^ = rfi- [(a + /3)<; ->!>“]+ . (35) 

We can express the Christoffel symbols in terms of the helicity-0 mode as 

rjr = - S [(a + m- mA . (36) 

which together with the relation 

[(a + )j)y - /St;;] s [(a +m- m] = a"". (37) 


implies that = 0 and we are left with 


6(j) 


-A {/?[(a + p)6; - /3t']}, 


(38) 


We see immediately that the equation of motion for the helicity-0 mode in the 
decoupling limit remains second order in derivatives even in the presence of the 
non-minimal couplings between the scalar held n and the ehective metric, after 
taking into account the equation of motion for the matter held. 
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5 Cosmological background 

In this section we would like to compute the field equations of the theory (1) on 
cosmological backgrounds. In a future work we shall push forward these preliminary 
studies and investigate the cosmological perturbations. Let us consider a FLRW 
background for the metric g and /, 

dSg = —N{tYdt^ + a{tydx‘^, and ds'j = —. (39) 

with a{t) being the scale factor. The effective metric then takes the form 

+ aefr(f)Mf^ , (40) 

with the effective lapse function and the scale factor as 


^efr(^) = (y.N{t) + (3nf(t) and aes(t) = aa(t) + (3af(t). (41) 


We can further define the Hubble parameters H = a/a, Hf = df/aj and H^s = 
des/des- Let us first compute the resulting effective energy density and pressure for 
the field tt 


P" 

p-K 


6ai HLtt + 


iL, 


N, 


eff . 
-TT 


eff 


9ao 


2 -"eff • 2 
71 


aoa. 




3 -f^eff • 3 
TT 


A® 


+ 1^#) 


0,2 


TT 


H. 


+ 4^if) 


2A3 


N: 


eff 


—ai(7r(3if/giVeff + 2iLeff) + 7l'(2iLeff — tt^) + 

JVeff iVeff 


(42) 


(43) 


in terms of which the Einstein equations for the g metric can be expressed as 

3 3 

2 2 

-2MU2H + 3H^N)+m^M^saNes^ = (44) 

Similarly the Einstein equations for the / metric are given by 

aj aj 

m^M^ViVeff^ = (45) 

Qj ^ f 
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On the other hand the equation of motion for the held vr on the FLRW background 
yields 

^ - HM) + j 

18(3 f 77 ^ / \ \ 

+ — ^3iVeff(-ffggiVeff + Hes) — 2ifeff-^eff j + 2-^^7r7rJ = aiReff(46) 

where we have dehned Res = 6{HeQ+2H^gNes). Here we have quoted the background 
equations of motion for our Lagrangian. They are quite involved but in principle 
solvable. The phase map is very reach due to the numerous degrees of freedom. In a 
follow-up work we would like to study the dynamical system of cosmological solutions 
by means of a phase map analysis and explore the critical points of the cosmological 
equations and compare with numerical solutions. It would be also interesting to see 
if a bouncing solution at early times could be possible. Furthermore, we would like to 
investigate the cosmological perturbations on top of this general FLRW background 
and study the stability conditions. 

6 Discussion 

In this work we have studied the theoretical implication of non-minimal derivative 
couplings of the effective composite metric, proposed in [35], to matter helds. Instead 
of considering the entire class of Horndeski scalar-tensor interactions, we focused 
on a very specihc subclass of non-minimal derivative interactions. The motivation 
for this comes from the fact, that these interactions are naturally generated by 
the covariantisation procedure of the decoupling limit of massive gravity. There, 
the specihc interactions between the helicity-O and hehcity-2 modes generate these 
specihc couplings between a scalar held and terms that depend on the Riemann 
tensor and Ricci tensor..etc. These curvature terms live in the usual space-time 
metric. Here we have promoted these derivate interactions to those living in the 
composite ehective metric. We hrst computed the covariant equations of motion. 
In the same way as in the original proxy theory, the Einstein and the double dual 
Riemann tensor are divergenceless, which then only generates second order equations 
of motion but in terms of the composite ehective metric. As next we performed a 
detailed Hamiltonian analysis of maximally symmetric space-times and considered 
three diherent scenarios. In the hrst scenario we assumed that the second metric is 
hxed and we saw immediately that the non-minimal derivative coupling makes the 
Hamiltonian to become non-linear in the lapse. In the second scenario we allowed the 
second metric to be dynamical as well but without including an additional kinetic 
term. In this case the Hamiltonian becomes automatically linear in the lapses. 
Finally, in the third scenario we included an additional kinetic term for the second 
metric and the consequence of this case was a Hamiltonian non-linear in the lapses. 
Thus, we concluded that the only consistent way of non-minimal derivative couplings 
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is through the second scenario with the dynamical / metric without an additional 
kinetic term. This is in a very close analogy to the studies in [54] in the context 
of allowed kinetic terms. Making the second metric dynamical enriches the phase 
space and gives more freedom. However, this has to be done without the inclusion 
of a second kinetic term. The conclusions of the hrst scenario with hxed / metric 
coincides with what was obtained in [46] for the Horndeski class. It would be very 
interesting to push forward their analysis of cosmological perturbations to the case 
of the second scenario with the dynamical / metric with no further kinetic term. 
We shall investigate this further in a future work. Furthermore, we investigated the 
non-minimal derivative couplings in the decoupling limit and saw that the equation 
of motion for the helicity-O mode remains second order in derivatives after taking 
into account the equation of motion for the scalar held. Moreover, we computed 
the held equation on FLRW backgrounds and we shall investigate the numerical 
solutions together with the cosmological perturbations in a future work. 
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